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Abstract. We derive the modulus of continuity of a class of processes called Canonical 
Embedded Branching Processes (CEBP), recently introduced by Decrouez and Jones 
[3], and we establish their monofractal character. CEBP provide a rich class of pro- 
cesses, including the Brownian motion as a particular case. The techniques developed 
in this study follow the steps of Barlow and Perkins on Brownian motion on a Sier- 
pinski gasket PQ, though complications arise here since CEBP are not Markovian in 
general. 
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1. INTRODUCTION 

The local fluctuations of a process X can be described using the local Holder ex- 
ponent hx{t), defined as 

hx(t) '■= liminf — ^— log sup \X(u) — X(t)\. 
e->o loge lu _ t{<€ 

When hx(t) is constant all along the sample path with probability 1, X is said to be 
monofractal. In contrast, there exist processes whose Holder exponent behaves errati- 
cally, whereby in any interval of positive length we find a range of different exponents. 
For such processes, it is in practice impossible to estimate hxif) for all t, due to the 
finite precision of the data. Instead, we use the multifractal or Hausdorff spectrum 
D x , a global description of its local fluctuations. D x {h) is defined as the Hausdorff 
dimension of the set of points with a given Holder exponent h. For monofractal pro- 
cesses, D x {h) degenerates to a single point at some h = H (so D X {H) = 1, and the 
convention is to set Dx(h) = — oo for h ^ H). When the spectrum is non trivial for a 
range of values of h, the process is said to be multifractal. 

Recently Decrouez and Jones [3] described a new class of processes, called Canon- 
ical Embedded Branching Process (CEBP) processes. CEBP are defined using the 
crossing tree, an ad-hoc space-time description of the process, and are such that the 
spatial component of their crossing tree is a Galton- Watson branching process. The 
contribution of this study is to obtain the modulus of continuity for CEBP processes, 
and hence show that they are monofractal. 
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Figure 1: A section of sample path and levels 3, 4 and 5 of its crossing tree. In the top 
frame we have joined the points T£ at each level, and in the bottom frame we have 
identified the fc-th level n crossing with the point (2 n ,T^_ 1 ) and linked each crossing 
to its subcrossings. 

2. THE CANONICAL EMBEDDED BRANCHING PROCESS (CEBP) 

Let X : M + — > R be a continuous process, with X(0) = 0. For nfZwe define level 
n passage times Tg by putting Tq = and 

T^ +1 = infO > T£ | X(t) G 2«Z, X(t) ± X(T£)}. 

The k-th level n (equivalently scale 2 n ) crossing C£ is the sample path from T^_ 1 to 
T" That is, CI = {(t,X(t)) | Tg_ x < t < T£}. 

When passing from a coarse scale to a finer one, we decompose each level n crossing 
into a sequence of level n — 1 crossings. To define the crossing tree, we associate nodes 
with crossings, and the children of a node are its subcrossings. The crossing tree is 
illustrated Figure [TJ where the level 3, 4 and 5 crossings of a given sample path are 
shown. 

In addition to indexing crossings be their level and position within each level, we 
will also use a tree indexing scheme. Let be the root of the tree, representing the 
first level crossing. The first generation of children (which are level —1 crossings, 
of size 1/2) are labelled by i, 1 < i < Zq, where Zq is the number of children of 0. 
The second generation (which are level —2 crossings, of size 1/4) are then labelled ij, 
1 < j < Zi, where Zj is the number of children of i. More generally, a node is an 
element of U = U m >oN m and a branch is a couple (u, uj) where u e U and j G N. The 
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length of a node i = %\ . . . % m is |i| = m, and the fc-th element is i[k] = i^. If |i| > m, 
i| m is the curtailment of i after m terms. Conventionally |0| = and i 1 = 0. A tree T 
is a set of nodes, that is a subset of U, such that 

• G T 

• If a node i belongs to the tree then every ancestor node i|fc, k < |i|, belongs to 
the tree 

• If u G T, then uj G T for j — 1, . . . , Z u and uj g" T for j > Z u , where Z u is the 
number of children of u. 

Let T m be the m-th generation of the tree, that is the set of nodes of length m. 
(These are level — m crossings, of size 2~ m .) Define T ; = {j G T | |j| > |i| and j||i| = i}. 
The boundary of the tree is given by <9Y = {i G M N | Vm > 0, i| m G T}. Let ipii) be the 
position of node i within generation |i|, so that crossing i is just C^Y\. The nodes to the 

left and right of i, corresponding to the crossings C^Yl^ and C^Y\ +1 , will be denoted 
i— and i+. In general when we have quantities associated with crossings we will use 
tree indexing and level/position indexing interchangeably. So Z\ = Z^YY T ; = T^Y), 
etc. 

Let a£ G {+, — } be the orientation of C%, + for up and — for down, and let be 
the vector given by the orientations of the subcrossings of C^. Let = T£ — T ( J i _ 1 be 
the duration of C%. Clearly, to reconstruct the process we only need a% and for all 
n and k. The a£ encode the spatial behaviour of the process, and the D% the temporal 
behaviour. Our definition of an EBP is concerned with the spatial component only. 

Definition 1. A continuous process X with X(0) =0 is called an Embedded Branching 
Process (EBP) process if for any fixed jiGZ, conditioned on the crossing orientations 

the random variables A% are all mutually independent, and A^ is conditionally 
independent of all A™ for m > n. In addition we require that {A% \ = i} are 
identically distributed, for i — +,—. 

That is, an EBP process is such that if we take any given crossing, then count 
the orientations of its subcrossings at successively finer scales, we get a (supercritical) 
two-type Galton- Watson process, where the types correspond to the orientations. 

Subcrossing orientations have a particular structure. A level n up crossing is from 
k2 n to (k + l)2 n , a down crossing is from k1 n to (k — l)2 n . The level n — 1 subcrossings 
that make up a level n parent crossing consist of excursions (up-down and down- 
up pairs) followed by a direct crossing (down-down or up- up pairs), whose direction 
depends on the parent crossing: if the parent crossing is up, then the subcrossings end 
up-up, otherwise, they end down-down. Let Z^ be the length of A%, that is, the number 
of subcrossings of C%. The number of up and down subcrossings will be written Z^ + 
and Z%~ respectively. Clearly, each of the Z£ — 2 first entries of A% come in pairs, each 
pair being up-down or down-up. The last two components are either the pair up-up or 
down-down, depending on Thus, given a£ = +, we must have Z£ + = \Z^ + 1 and 
Z^~ — -xZ 1 ^ — 1, and conversely given a% — — . 

Let A be the space of possible orientations. That is, a G A consists of some 

number of pairs, H — or — h, then a single pair ++ or . Given an EBP process, 

for the offspring type distributions we write p\(a) = P(v4£ = a | a% — +) and p^(a) = 
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P(A£ = a | al = -)> for a G A Let //+ = E(Z™ | a£ = +), /i" = E(Z™ | a£ = -) and 
/i = + then the mean offspring matrix is given by 

(^ n+ K = -) (^ri< = -w U^- 1 K + i; 

To proceed we need to make some assumptions about pj. 

Assumption 1. fi~ > 2 and E(Z™ log | a£ = j) < oo /or i, j = ± 

The first of these assumptions ensures that M is strictly positive with dominant 
eigenvalue /i > 2, and corresponding left eigenvector (|,|). The corresponding right 
eigenvector is — 2)/(/i — 2), — 2)/(/i — 2)) T . The second assumption is the 
usual condition for the normed limit of a supercritical Galton- Watson process to be 
non-trivial. 

Theorem 1. (f^j Theorems 2.1 and 2.2) For any offspring orientation distributions 
Pa satisfying Assumption [21 there exists a corresponding continuous EBP process X 
defined on R + . Moreover, we can choose X so that for each n £ Z, conditioned on 
the crossing orientations a%, the crossing durations D% are all mutually independent, 
and D% is conditionally independent of all A™ for m > n. Also, E(D^ | a% = ±) = 
/^ n (/i ± — 2)/(/i — 2), and the distribution of ji~ n D^ depends only on a%. 

Up to finite dimensional distributions, X is the unique such EBP with offspring 
orientation distributions . That is, for any other EBP process Y with offspring ori- 
entation distributions p^ and crossing durations as above, we have (X(ti), . . . , X(tf.)) = 
(Y(ti), . . . , Y(tk)) for any < t\ < t% < ■ ■ ■ < tf.. Accordingly we call X the Canonical 
EBP (CEBP) process with these offspring distributions. 

The two-type Galton-Watson process defined by the orientation distributions p^ is 
supercritical, and when scaled by its mean converges to (~, \)W l , where i = ± is the 
type of the first individual. For the Canonical EBP we have that the crossing duration 

D% = fi n W\ where i = a%. 

We also observe that X is discrete scale-invariant: let H = log 2 / log \i, then for all 
c e £ Z} ; 

X(t) f = c~ H X(ct), (1) 

where = denotes equality for finite dimensional distributions. H = log /i/ log 2 is 
known as the Hurst index. 

To simplify our lives somewhat, we will from here on restrict ourselves to EBP 
for which the embedded Galton-Watson process is single type. That is, the number 
of subcrossings does not depend on the orientation of the parent crossing. We will 
also suppose that, given the number of subcrossings, the excursions are independent 
and equally likely to be up-down and down-up. That is, we will make the following 
assumption. 

Assumption 2. The distribution of is independent ofa%, and 

P + a(- •• + +)= Pa( ) = 2-*P(Z? = 2{z + 1)), 

where ■ ■ ■ represents a combination of z pairs, each either H — or — h 




4 



Write Z for a generic r.v. with the same distribution as the Z£ 
When Assumptions [Hand [2] hold, the CEBP is completely specified by the subcross- 
ing family size distribution Z . Starting with any single crossing, if we just count the 
number of subcrossings at successively finer scales, then we get a single-type Galton- 
Watson process with offspring distribution Z . Clearly /i + = \x~ = \x = KZ. Let W be 
the limit of the (supercritical) Galton- Watson process with offspring distribution Z, 
normed by its mean, then KW = 1 and for any given n the crossing durations D% are 
all distributed as fi n W. 

3. THE MODULUS OF CONTINUITY OF A CEBP 

The goal of this section is to establish that CEBP are monofractal processes with 
Holder exponent H = log 2/ log /x £ (0,1). Our approach is based on that used by 
Barlow & Perkins pQ to obtain the modulus of continuity of Brownian motion on 
the Sierpinski gasket, though complications arise because CEBP are not in general 
Markovian. 

The basic idea is to use bounds on the crossing durations D% to control how fast the 
process can move away from a given point. We will take Assumptions [1] and [2] to hold 
throughout. Let X be the CEBP determined by the subcrossing number distribution 
Z. 

We have the following from Biggins & Bingham [2J. 

Lemma 1. Suppose that Assumptions U\ and hold. There exists strictly positive 
constants c\, c-i, and C3 such that for all x > 

exp {-cix~ HI{1 - H) ) < F(W < x) < c 2 exp (-c 3 x~ H/(1 ~ H) ) . 

Let Tj}(s) be the k + 1-st level n crossing time greater than or equal to s, for k > 0. 
So if s is a level n crossing time then Tq(s) = 0. The previous lemma gave us a 
bound on the duration of a crossing. The next lemma gives a lower bound on the 
time remaining in the current crossing. To establish this result we will need to make a 
further modest restriction to the class of CEBP we consider. 

Assumption 3. We assume that the subcrossing number distribution Z is such that 
there exists a £ such that for all y 

Z + C >st Z - y I Z > y. 

Here > st denotes stochastic domination. That is, for all y and z, 

F{Z - y > z \ Z > y) <F(Z + ( > z). 

This condition clearly holds for Z bounded, and for Z that are NBU (New Better 
than Used, in which case £ = 0). Examples of NBU distributions include the negative 
binomial with shape > 1 and the Poisson. 

Let {J r s } s >o be the filtration generated by X. 

Lemma 2. Suppose that Assumptions [21 [H and [3| hold, then there exists constants 
C4, c 5 > such that for all x > and n £ Z, 

P(T ( ; t (s) ~s<x\F s )> c 4 exp (-c 5 (^- n x)- H/{1 - H) ) . 
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Figure 2: Notations in Lemma HJ The current level n crossing at time s is in bold. The 
first level n crossing time greater than s, denoted Tq(s), corresponds to T%, and the 
first level (n + 1) crossing time greater than s is Tq +1 (s) = T£. There are Z n+1 (s) = 3 
level n crossings from Tq(s) to Tq +1 (s), and y — 1 level n crossing from the current 
level (n+1) crossing that has already happened at time s, including the current crossing 
in bold. 



Proof. Note first that 

Z n + 1 {s) 

T n+1 (s) = T n (s)+ 

i=l 

where Z n+l (s) > is the number of level n crossings from T n (s) to T n+1 (s), and the 
W(i) are independent and distributed as W. If Z n+1 (s) is not zero then, conditioned 
on J-* s , it will be distributed as Z — y | Z > y, where Z has the subcrossing number 
distribution, and y > is the number of level n crossings from the current level n+1 
crossing which have already happened by time Tq(s), including the current level n 
crossing. Notations are illustrated in Figure [2J 

Thus, from our assumption on Z, conditioning on T s we have 

z+C 

T- + \s)-T-(s) < st ^T//W(*) 

i=i 

C 

1=1 

As n 4- — oo we have Tq(s) 4 s (this follows directly from [1] Theorem 1), whence 

n-l c 

TZ(s)-s< stf i n W(0,n)+ zZ^ W ^ k ^ 

k=— oo i=0 
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where the W(i, k) are i.i.d. with distribution W. Thus, for any 9 > we have 

n— 1 , 

E(e~ 9iT ^-^\T s )>Ee^ nw J] (le"^ V . (2) 

fc=— oo 

Using the lower bound in Lemma [1] for the left tail of W, it follows that for any x > 

> exp (-Ox - Cl (xfi~ k y H/il ~ H] ) . 

For x = (c 1 /9) 1 - H fi kH we get Ee"^ few/ > exp (-c 2 (9j2 k ) H ) . Plugging this into © 
yields the bound 

E(e-^W-) | F s ) > exp (-c 3 (^)^) . 

Applying Lemma 4.1 in [T] with 9 = (c4/x) 1 ^ 1_H ^2^ 2+fc ^ ( ^ 1 ~ H \ and readjusting the 
constants, gives the result. □ 

Lemma 3. Suppose that Assumptions [21 and hold, then there exist constants 
Cq, . . . , Cg > st/cn that for all A > and any s, t > 0, 

c 6 exp(-c 7 (A 1/ ^/t) H/(1 -^) <F(\X(s + t)-X(s)\ > A | F s ) 




Proof. We start with the last inequality. Define n E 7h by 2 n < A < 2™ +1 . If the 
maximum variation of X in the interval [0, t] is at least A, then necessarily there exists 
a k such that s < T£Z\ < T£~ l < s + t. Thus, using Lemma [fl 

P ( sup \X{s + u) - X(s)\ > A | JF 8 ] < Ppr 1 < <) 

\0<u<t / 

<c 2 exp(- C3 (t/^ 1 )^/( 1 ^) (3) 

Re-expressing the last inequality in terms of A and adjusting the constants yields the 
desired upper bound. 

We now turn to the first inequality of the lemma. Our proof is based on Theorem 
4.3 in prj, though more work is required because our process is not Markov. This time 
let n e Z be such that 2 n ~ 2 < A < 2 n ~\ As before let T fc n (s) be the k + 1-st level n 
crossing time of the process after time s, and also let T!! 1 (s) be the first level n crossing 
time strictly before s. 

Consider the possible level n movements of the process up to time Tq(s). We take 
cases depending on the orientations of the two level n crossings leading up to Tq(s). 

(a) (b) ++ (c) +- (d) -+ 

Define tt = P(Z > 2) (> by Assumption Q} . 

In case (a) the next two level n crossings will have orientation — + with probability 
7r/2. Let a" x (s), Oq(s), a™(s), a^s), be respectively the orientations of the two cross- 
ings up to time Tq(s) and the two crossings after time Tq(s). Also let T™ _1 (T 1 ri (s)) be 
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Figure 3: The points considered in case (a) of the proof of Lemma [3j 

the next level n — 1 crossing time after T™(s). Let Dq(s) = Tq(s) — T™ 1 (s) = fi n W, 

D?(s) = Tf(s) - T n (s) = fi n W and D? _1 (77 l (s)) = Tf-^T^s)) - T?(s) = . 
Note that they are independent and independent of the a£(s). We have 

F(\X(s + t)-X(s)\ > A | T s ,a\(s) = -a%(s) = -) 

> F(\X(s + t)- X(s)\ > 2"" 1 | T s , al^s) = - etf(s) = -) 

> n/2¥(\X(s + t) -X(s)\ > 2 71 - 1 \ F s , 

a\(s) = -,a n (s) = -,aUs) = -,a n 2 (s) = +) 

> n/2F(T?{s) < 8 + t, T?-\T?(s)) >s + t\F s , 

al^s) = -,a%(s) = -,a?(s) = - a£(s) = +) 

> tt/2P(T "(s) < s + t/2, D?(s) < t/2, Dr\T?(s)) >t\JR s , 

a n _ x {s) = -,a$(s) = = -,a£(s) = +) 

= tt/2P(T "(s) < 8 + t/2 | j;, a^s) = -, a»(s) = -) 
F(D?(s) < t/2)P(^" 1 (T 1 "(s)) > t) 

Thus from Lemmas [T] and [2] we have 

P(|X(s + t) -X(s)\ > A | J- aj a^(5) = = -) 

> c a exp (-c 2 (A 1/ - ff /t)' ff/(1_H) ) (1 - c 3 exp (-^(A 1 /*/*)^ 1 "^)) 

Choose if large enough such that the last term is at least 1/2 when A 1 ^ jt > K. Thus, 
since the LHS is decreasing in A, we can find C5 such that for all t G [0, 1] and A > 0, 
the LHS is bounded below by c 5 exp (-^(A 1 / 2 */*)^ 1-2 ^). 
Cases (b) is analogous to case (a). 

In case (c) we distinguish two further possibilities: (cl) the next two level n crossings 
form an excursion (either — h or H — ); and (c2) the next two level n crossings form 
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a direct crossing (either or ++). In case (cl) with probability 1/2 the excursion 

will be — h, in which case we can proceed as in case (a) to get a bound of the same 

form. In case (c2) if the direct crossing is then the approach of case (a) again 

suffices, however if the direct crossing is ++ then we need to modify the argument a 
little. In this case we wish to bound ¥(\X(s + t) —X(s)\ > A | J r s ,a™ 1 (s) = +, cKq (s) = 
— , a" (s) = +, «2 (s) = +). With probability 7r/2 the next pair of level n crossings are 
the excursion H — . Let T™~ l {T£(s)) be the next level n — 1 crossing time after T£(s), 
then we get 

¥{\X{s + t)-X{s)\ > A | F^a^s) = = - = = +) 

> tt/2P(|X(s + 1) - X(s)\ > 2 n ^ | JF S , 

al^s) = +,a% >s = -,a?(s) = +, a£(s) = = = -) 

> tt/2P(T 3 "( S ) < s + ^Tr 1 ^)) > s + t | T s , 

al^s) = = = = = +, = -) 

> vr/2P(T "( S ) < s + t/4, D?(s) < t/4, D n 2 {s) < t/4, D%(s) < t/4, D?" 1 > * I T> 

al^s) = +,a£, = -,a?(s) = +,a£(s) = = +,<(s) = -) 

= 7r/2P(T "(s) < s + t/2 I T s , a^s) = +, = -) 

P(D"(s) < i/4)P(D£(s) < */4)P(I£(s) < t/4)P(^- 1 (T 3 n (s)) > t). 

This can be bounded below in the same way as in case (a). 
Case (d) is analogous to case (c). 

Finally, for general t > 0, let m be such that fi~ m t < 1. Then, noting that 
{2~ rn \) 1 ' H /(tfi~ m ) = X^ H /t > K, by the discrete scaling of X, 

F(\X(s + t) -X(s)\ > A | T s ) = P(|X(/i~ m (s + 1)) - X(/i" m s)| > 2~ m A | J>-™ s ) 

> c 5 exp (-c 6 ((2- m A) 1/ ^/(/i- m t)) H/(1 - H) ) 
= c 5 exp(- C6 (A 1 ^/t)^/(^)), 

which concludes the proof of the lemma. □ 

Theorem 2. Suppose that Assumptions U\ IE and\^ hold. Let = 5 H | log <5| 1— ^ 7 

then there exist constants cio, ch > sitc/i i/itrf 

Cio < hmmt sup 



5 -> s,*6[0,l],|t-s|<5 ^ff(* — s ) 

\X(t)-X(s)\ 
< hmsup sup — - — -; r < C\\. 

5^0 s,te[0,l],\t-s\<6 h H (t — s) 

Proof. Consider first the lower bound. Fix c\ > 0, then for any I > and m = 
0,1,..., 2' - 1, put 

A m j = {\X((m + l)2~ l ) -X{m2- l )\ > c l l 1 - H 2~ lH } . 
By Lemma [3] we have F(A m j \ J- m 2-i) > C2e~ C3 ', where C3 oc c^ 1 . By repeatedly 
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conditioning we have 

(2 ! -l \ 2 l -l 2 l -l 

n A ki = n F t A ki i -^-o = n ^ - p ^™.< i ^-o) 
m=0 / m=0 m=0 



< c 4 exp (-c 2 e (log2 - C3)/ ) 



We can choose Ci so that log 2 — c 3 > 0, in which case the RHS above tends to as 
/ — > oo, and we have 

P (\X(t + 2~ l ) - X(t)\ < c 5 h H (2~ l ), VZ > 0, t G [0, 1 - 2"']) = 0, 

which establishes the lower bound. 

For the upper bound we proceed in a similar manner, though we can no longer 
just consider points on the lattice 2~ l 7L. For / > and m = 0, . . . , 2 l — 1, let = 
[m2~\ (m + l)2~ l ), and define 

<D m ,, = sup \X(t) - X(m2- l )\ 
B m ,i = {<S> m ,i>c 1 l 1 - H 2- lH } 

From our estimate in Lemma [3] we have F(B m j \ J 7 ^- 1 ) < C2e _C3 ', where C3 oc cj^ 1 . 
Thus, repeatedly conditioning on T m 2-i for m = 2 l — 1, . . . , 0, we have 



F(B m>l for some < m < 2 l ) = 1 - ¥(B^ for all < m < 2 l ) 

< l-(l-c 2 e" 



-c 3 l\2 l 



c 2 e" C3 '2' N 



2 l 

< 1 -exp (-c 4 e- (C3 - log2)/ ) 

< c 4 e- (c3 - log2)i . 

Here we have chosen c\ so that C3 — log 2 > 0. 

Applying the Borel-Cantelli lemma, we see that there exists an L such that with 
probability 1 

$ m j < c 1 l 1 - H 2- lH for all / > L and < m < 2 l . 

Now let s G I m ,i and suppose that t is such that s < t and \s — t\ < 2~ l . Then 
t G J m j U Im+1.1 and we have, with probability 1, 

\X(t)-X(s)\ < \X(t)-X((m + l)2- l )\ + \X((m + l)2- l )-X(m2- l )\ 

+ \X(m2- 1 ) -X(s)\ 
< 3c 1 / 1 ^2~^. 

If we take 2~^ +1 - ) < 5 < 2~ l , then we have, with probability 1, 

sup \X(s) - X(t)\ < c 5 l 1 - H 2- lH < c G h H (5), 

s,t£[0,l],\s-t\<8 

as required. □ 
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Remark 1. In the special case where the CEBP reduces to a Brownian motion, the 
existence of the limit in Theorem^ follows from Levy's modulus of continuity theorem. 

Corollary 1. Suppose that Assumptions^ [H and\^hold, then the CEBP is a monofrac- 
tal, in that P-a.s. the Holder exponent h(t) = H for all t G [0, 1]. 
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